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We study Casimir interactions between cylinders in thermal non-equilibrium, where the objects 
as well as the environment are held at different temperatures. We provide the general formula for 
the force, in a one reflection approximation, for cylinders of arbitrary radii and optical properties. 
As is the case for equilibrium, we find that the force for optically diluted cylinders can be obtained 
by appropriate summation of the corresponding result for spheres. We find that the non-equilibrium 
forces are generally larger than their equilibrium counterpart at separations greater than the thermal 
wavelength. They may also exhibit oscillations as function of separation, leading to stable points of 
zero net force. These effects are particularly pronounced for thin conducting cylinders (e.g. 40nm 
diameter nano-wires of tungsten) due to their large emissivity. 
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I. INTRODUCTION & SUMMARY 

There has been considerable recent activity in fabrica- 
tion of micro and nano electromechanical devices which 
actuate forces in response to currents. At sub-micron 
scales classical electromechanical forces are supplanted 
by forces due to quantum and thermal fluctuations pQ. 
In equilibrium these are manifested as van der Waals and 
Casimir forces [5] , while the passage of currents may well 
lead to temperature imbalances resulting in thermal ra- 
diation and non-equilibrium forces. Fluctuation-induced 
forces are typically non additive, and (especially for con- 
ductors) with sensitive dependence on shape and mate- 
rial properties. Nanowires and nanotubes provide a sim- 
ple example of an extended shape, with several possible 
applications [3HS]- Moreover, cylindrical geometries are 
amenable to high precision experiments [7Ml2j. provid- 
ing good contrast to the more widely studied spherical 
geometry |13j . Cylinders can also be easier fixed posi- 
tionally and held at different temperatures for possible 
experimental study of non-equilibrium effects [14] . 

In his seminal paper |15| Casimir calculated the force 
between two perfectly conducting parallel plates, due 
to the quantum zero point fluctuations of the electro- 
magnetic (EM) field in the intervening vacuum. Lif- 
shitz [16 later extended Casimir's computations to the 
case of real dielectric materials and to finite tempera- 
tures. A key step in Lifshitz's approach is to include 
thermal and quantum fluctuations of currents (sources) 
in the dielectric, following the formalism of fluctuational 
electrodynamics pioneered by Rytov j!7j . Typically at 
small separations zero-point fluctuations shape the force, 
whereas at separations large compared to the thermal 
wavelength At, thermal effects dominate [THl HH1 HI] • fty- 
tov's formalism is also appropriate to out of equilibrium 
steady states in which each object is held at a different 
temperature. There is extensive literature on the topic 
of non-equilibrium interactions between two atoms (or 
molecules) 20 20 . Recently out of equilibrium Casimir 
forces have been considered in several systems: parallel 



plates [27H29] . modulated plates [30], plate and an atom 
in different setups [551 I3"TH3~4"] . two spheres and sphere 
and a plate [18) . Formalisms for arbitrary objects were 
presented in Refs. [2H1ES]- A common feature for setups 
involving compact objects is the need to account for the 
contribution of the environment to the force, which de- 
pends on a possibly different ambient temperature. 

In this paper, we consider non-equilibrium forces be- 
tween two parallel cylinders characterized by an arbitrary 
(including axially anisotropic) dielectric function. Us- 
ing the Rytov formalism general expressions are obtained 
for forces between wires maintained at different temper- 
atures from each other and the environment. These ex- 
pressions are then analytically and numerically studied in 
a number of cases. With the general reader in mind, in 
the remainder of this section we provide both an outline 
of the paper as well as a summary of its main results. 
The interested reader can then proceed to the detailed 
derivations that follow in subsequent sections. 

Analytical results are presented in Sec.[nJ starting with 
a brief introduction to the Rytov formalism in |II A[ The 
fluctuating EM field in the space between the cylinders is 
sourced by fluctuating currents in the cylinder, with an 
additional contribution from the surrounding. The cur- 
rent correlations in each source can be related by a fluc- 
tuation dissipation theorem to the corresponding tem- 
perature |17j : the stress tensor (from which forces are 
computed) is then related (via Green's functions) to the 
source fluctuations. The force on each cylinder now has 
two contributions: an interaction force sourced by the 
other cylinder, and a self-force due to modification of the 
EM fields sourced by itself through influences of the other 
cylinder. These two contributions are computed respec- 
tively in sections |IIB| and |II C[ within a so-called one 
reflection approximation in which multiple scatterings of 
the EM field sourced by either cylinder are ignored. The 
results are further simplified for thin cylinders, with radii 
R\ and i?2 small compared to thermal wavelength, rele- 
vant skin-depth and separation d. In this limit, the non- 
equilibrium force per unit length is proportional to R\R\, 
and decays at large separation as It is well known 
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FIG. 1: (color online). Total force on cylinder 1 per unit 
length in a system of two SiC cylinders with equal radii 
R = O.lpim at separation d in a a) cold (0 K) b) warm 
(300K) environment. Dashed lines indicate repulsion. Points 
of change from repulsive to attractive force with increasing d 
correspond to stable points of zero force. 



FIG. 2: (color online). Total force on cylinder 1 per unit 
length in a system of two tungsten cylinders with equal radii 
R = 0.02/im at separation d in a a) cold (0 K) b) warm 
(2400K) environment. Dashed lines indicate repulsion. Hor- 
izontal line Fg/L denotes the weight of the corresponding 
tungsten cylinder per unit length. 



that in the optically dilute limit, with dielectric response 
c — > 1, equilibrium Casimir forces become pair- wise addi- 
tive, obtained by summing over contributions of pairs of 
polarizable volume elements. In Sec. |HD| we show that 
a corresponding result holds out of equilibrium, namely 
that the force on cylinders of sufficiently diluted media 
can be obtained by pairwise addition of corresponding 
forces for chains of small spheres [TB] by treating a thin 
cylinder as a chain of small spheres. 

As discussed in Sec. precise calculations of the 
force can be performed numerically for different materi- 
als, starting from optical data in the form of a (frequency 
dependent) dielectric response s(u>). In particular, we re- 
port results for silicon-carbide (insulator) and tungsten 
(metal): Figure [l] depicts forces for two SiC nanowires of 
radii R\ = i?2 = 0.1/im; each wire or the environment 
is at OK or 300K, for a total of 8 possible combinations. 
The top panel correspond to the four cases where the en- 
vironment is at OK, the bottom to where it is at 300K 
ambient temperature. To depict results on a logarith- 
mic scale, attractive forces are indicated by solid lines 



and repulsive forces by dashed lines. The equilibrium 
force (whether at OK or 300K) is attractive at all separa- 
tions d] non-equilibrium forces deviate strongly from the 
equilibrium force at distances larger than the thermal 
wavelength \x(300K) w 7.6/xm. The most significant 
deviation is for the force on a cold cylinder due to a hot 
cylinder at T env = OK (top panel) which becomes repul- 
sive due to radiation pressure, falling off as 1/d at large 
separations. Note that for the same setup the force ex- 
perienced by the hot cylinder is not equal and opposite, 
instead undergoing oscillations at a wavelength related to 
the resonance of SiC. Oscillations are in fact a common 
feature in these curves, which thus go through several 
points of zero net force, alternating between stable and 
unstable mechanical equilibria. 

Conducting cylinders strongly radiate for radii of 
the order of the relevant skin-depth; hot conducting 
nanowires are hence good candidates for experimen- 
tal observation of thermal nonequilibrium effects. Fig- 
ure [2] depicts the corresponding force curves for tungsten 
nanowires of R\ = Ra = 20nm. The higher tempera- 
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ture considered in this case is 2400K; the large temper- 
ature and small radius are chosen for emphasis. Indeed 
in this case both the equilibrium and non-equilibrium 
forces are six orders of magnitude larger than for SiC 
(non-equilibrium effects are roughly three orders of mag- 
nitude larger for tungsten at room temperature compared 
to SiC, and additionally enhanced by increasing the tem- 
peratures of the tubes). The overall trends (repulsive 
or attractive forces) are similar to SiC, but there are 
fewer oscillations. Interestingly, for conductors, the non- 
equilibrium correction to the Casimir forces can be two 
orders larger than the weight of the cylinder, thereby be- 
ing in the observable range in principle. We also give a 
rough estimate of Ampere's force between the nanowires, 
based on the currents necessary to keep their temper- 
atures constant, which is much smaller than the non- 
equilibrium correction, rendering conducting nanowires 
potential candidates for experimental studies of non- 
equilibrium Casimir forces. 1 

The above results (as well as input dielectric func- 
tions and other details) for SiC and W are discussed in 
Sees. Ill A and |III B] Additional computational details, 
such as the expansion of scattering amplitudes at small 
cylinder radius, are relegated to Appendices. 



II. NON-EQUILIBRIUM CASIMIR FORCES 
BETWEEN CYLINDERS 



A. Formalism 

We apply the formalism introduced in Ref. [35] and 
further discussed in Ref. [TB] , to a system of two parallel 
cylinders of radii Rj (j = 1, 2) and axis-to-axis separa- 
tion d in vacuum. Each object is held at a constant and 
homogenous temperature {T^}, and embedded in an en- 
vironment at temperature T env . The cylinders are char- 
acterized by electric and magnetic responses Ej and fij. 
For anisotropic materials Sj can be a tensor; we special- 
ize to objects with azimuthal and translational (in the 
direction of the cylindrical axis) symmetries such that 
the problem can be decomposed into independent sec- 
tors indexed by n and k z (see below) . Such a description 
is valid for isotropic as well as for uniaxial materials (with 
optical axis parallel to the cylinder axis); such a model 
was recently proposed as a description of multi-walled 
carbon nanotubes 37J. 

In a non-equilibrium steady state, the net force on the 
system may be non-zero since there is momentum carried 
away by the field radiated to the environment. Thus we 



1 We note however that the validity of a continuous dielectric func- 
tion when the size of the object is comparable to the mean free 
path of the charges (typically a few nanometers for conductors) 
is an issue to be addressed in future studies. For example, plas- 
mon resonances can occur in this case 36_, which might change 
the non-equilibrium phenomenology in interesting ways. 



need to consider the force acting on each cylinder sepa- 
rately. We denote the total force acting on the cylinder 
1 by F' 1 ^ , whereas F^ 2 ) can be found directly from the 
expression for F^ by interchanging indices 1 and 2 and 
changing its sign. As shown in Refs. [T51 [33] . 



(1) 



j=l,2 



where F^- eq (T env ) is the Casimir force between the 
cylinders for the case of global equilibrium at temper- 
ature T env , containing the contribution from zero point 
fluctuations. This force is not discussed in the present 
paper and treated as known (it can be computed using 
scattering results for cylinders as in Refs. [Sl|3Hj; more 
general numerical schemes exist for calculating Casimir 
forces between 3D objects of arbitrary shapes and dielec- 
tric properties 39J). The difference of F^ (T env , T x , T 2 ) 
from F^' eq (T env ) is due to the deviations of the cylin- 
der temperatures T\ and T 2 from T env . F^ (T) is the 
force acting on cylinder 1 due to the sources in cylin- 
ders j = 1,2 at temperature T. Although dealing with 
three sources (including the environment), we have thus 
only to evaluate two terms: F^(T) and Fy(T). Due to 
symmetry of the system, the net force is parallel to the 
axis-normal connecting the two axes, and we denote its 
magnitude by F X (1) (T) and F 2 (1) (T). 

Computation of forces proceeds through integrating 
appropriate components of the EM stress-energy ten- 
sor around a surface enclosing one (or both) cylinders. 
In the Rytov formalism [17', correlations of the electric 
field E at frequency oj are sourced by sources in the two 
cylinders. The contribution of each cylinder can be com- 
puted separately [37] , and for a pair of points at r and r' 
outside the cylinder, the symmetrized correlator has the 
following spectral density, 

/OO J 

C j (T j , U ) = -a(T j ,u>) Yl E / 
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dk z 
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where M ni £ z and N B] ), ( are the two polarized outgoing 
cylindrical waves, see App. [A] These waves are indexed 
by the multipole order n and k z , the component of the 
wavenumber k — uj/c along the cylinder axis. The overall 
strength of correlations is set by 

a(T,cj) = uj2h ^ 2 (exp[^/fc B T] - l)' 1 , (3) 

which describes the occupation number of all oscillators 
of frequency w; c is the speed of light, and H is Planck's 
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constant. We have also introduced the abbreviation 
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describing the contribution of propagating and evanes- 
cent waves. Note that only propagating waves contribute 
to the heat emitted by a single cylinder, whereas evanes- 
cent waves also contribute to the interactions of two 
cylinders. T^ k ^ is the T-matrix element of cylinder j, 
which relates the amplitude of a scattered wave of polar- 
ization P in response to an incoming wave of unit ampli- 
tude and polarization P' . The explicit form of T-matrix 
elements for isotropic or uniaxial materials can be found 
in Sec. Ill of Ref. [57] (see also [3D] for the isotropic case). 
In general, the T operator of a cylindrical object is sym- 
metric and diagonal in n and k z , but couples different 
polarizations. 



B. Interaction Force 

In this subsection we compute the force exerted on 
cylinder 1 by the field produced by cylinder 2, i.e., the 



interaction force F2 (T). To this end, we scatter the 
field in Eq. ^ (for j = 2) at cylinder 1 and then com- 
pute the force by integration of the Maxwell stress tensor 
on a surface enclosing 1. Since multiple reflections on the 
cylinders are ignored this is a one reflection approxima- 
tion [18], valid for large d/R, which results in a force per 
unit length of 
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(5) 



where q = yjk 2 — k 2 is the wavenumber component per- 
pendicular to the cylindrical axis. Note that negative 
and positive signs of the numerical result correspond to 
attraction and repulsion respectively (again, the force is 
parallel to the axis-to-axis separation). The functions / 
and / in Eq. ^ take the rather lengthy forms 
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where stands for the Hankel function of first kind 
of order n. The interaction force in Eq. |5| consists of 
two distinct terms: Due to propagating (|fc7[ < uj/c) and 
evanescent waves (\k z \ > w/c), whose properties will be 
discussed below. We note that for nanowires with thick- 
ness in the range of a few tens of nanometers, Eq. ([5| 
gives accurate results down to separations of well below 
1 micron. 

Thin cylinders ({Rj} <C d, {Sj}, At): The general ex- 
pression in Eq. ^ is complicated as it involves an in- 
finite series; more insightful expressions are obtained in 
the asymptotic limit of thin cylinders. In this paragraph 
we analytically study the limit {Rj} <C d, {Sj}, At, where 



c/Im[y/Ej]u! is the skin-depth of tube j. The ad- 



ditional limits {Rj} <C {5j},\t allow restriction to a 
finite number of partial waves whose T-matrix elements 
are proportional to R 2 (see App. [Bland note that we take 



li = l). Thus all terms with |n| or \m\ larger than 1, and 
the terms with products of two T-matrix elements, can 
be neglected. The corresponding force can be found ana- 
lytically in the two limits concerning the ratio of thermal 
wavelength and separation. For d <C At 2 , we have 



lim 
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L 



= h 
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H 1 U 2 



c 2 d 4 



(7) 



where the auxiliary functions g™ and g™, that depend 
only on the dielectric functions, are given in Eqs. (C.l) 
and ( C.2 ), respectively. We note that the force in Eq. (|7) 
is in most cases attractive and dominated by the evanes- 
cent part. (It can nevertheless be made repulsive for 
certain dielectric functions, as discussed in Ref. [UJ). At 
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large separations, i.e., d ^> At 2 , the leading term of the 
interaction force is, 



lim 
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where g\ n is given in Eq. (C.3|. This limit always yields 
a repulsive force which is due to momentum transfer by 
propagating waves. Evanescent waves do not contribute 
in this order as they decay too quickly with separation. 

Further insight can be gained by additionally requir- 
ing the temperature to be so low that one can expand the 
dielectric function, i.e., At 2 3> Aq, , where Ao^ is the wave- 
length of the lowest resonance of the dielectric response 
of cylinder j. According to Ref. [42] . 

lim Ej (u) = £o 3 + i^^- + 0{w 2 ), (9) 

Ao 3 <CAt 2 C 

with real valued static dielectric constant So j and inelas- 
tic collision length Aj„. . With this form of £j(uS), the 
corresponding leading behaviors of the force are 
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Note that the dependence of temperature (via powers of 
At) is quite different in these limits. The correspond- 
ing auxiliary function s f™ , /|™ and f[ n are given by 
Eqs. Jail), ((051) and dcB, respectively. 



C. Self- force 

We next compute the self force F^\ acting on cylin- 
der 1 due to the field emitted by cylinder 1 and reflected 
from cylinder 2, again restricting to d 3> {Rj}, where the 
one reflection approximation is valid. The origin of this 
force is cylinder 2 acting as a reflector of EM waves emit- 
ted by cylinder 1. The interference of the emitted and 
reflected waves is expected to give rise to the observed 
oscillatory behavior of the self-force as a function of sep- 
aration. Computationally it is easier to first evaluate the 



force F 



(1+2) 



(2) i 

F{ +F1 1 ' acting on both cylinders due to 

(2) 

sources on cylinder 1, and to then subtract F^ . Again, 
it is a property of non-equilibrium that the net force on 



(2) 

the system is not zero. Note that F^ is obtained from 
Eq. ([5 ) by interchanging indices 1 and 2 and changing its 
overall sign. 

After tedious algebra we arrive at the following expres- 



sion for the force F x 
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s i,n,m,k z = 2 Re A 1 £,k z I m 



H 



(i) 



l (qd)J n -,n-i(qd)T[^ kz 



-J n _ m (gd)i?™^ n _ 1 (gd)T; 



ipp'* 
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and J n denotes the Bessel function of order n. As ex- 
pected, the force on the system is solely due to propagat- 
ing waves, since evanescent waves do not carry momen- 
tum to the environment. 

Thin cylinders ({Rj} <C d, {Sj}, At): We gain further 

insight by examining i^ 1+2 ** in the asymptotic limit of 
thin dielectric cylinders. For d <C A^ , 



F (i+2) 

lim -± = OfcT 1 ) 

d,{ < 5 j },AT 1 >{ii J -} L 



(14) 



i.e. diverging weakly at small d. Compared to this, the 
force F^ in Eq. Q is proportional to d~ 6 , and we can 
neglect F± 1+2 ' in this regime, so that for d <C Ati 



lim = 



lim F> 

d,{5j},X Tl »{Rj} 



(2) 



We thus observe that at close separations these contribu- 
tions to the forces on the cylinders arc equal and opposite, 
again because momentum transfer to the environment 
can be ignored. 

In the opposite limit of large distances, d 3> A^ , one 
can verify that 



,(i) 



lim 



rf,{5 J },A Tl >{fi; 3 } L 



0(d~^ 2 ), 



(15) 



where the prefactor is a lengthy algebraic function ad- 
ditionally depending on d in an oscillatory manner via 
exp (2iujd/c). Thus, at large separations the self-force is 
asymptotically less relevant compared to F% (following 
d : in Kq. flii. 



D. Cylinders from spheres in the dilute limit 

A rarified material can be regarded as a collection of 
independent molecules; forces between two such materi- 
als may then be obtained by pairwise summation of forces 
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FIG. 3: Casimir forces between objects can also be obtained 
by pairwise summation of the forces between all the volume 
elements. 



between their molecules. Indeed, equilibrium Casimir 
forces in the optically dilute limit (sj(oo) — > 1) can be cal- 
culated by integrating the interactions between volume 
elements. In this subsection we show that even the non- 
equilibrium force between two cylinders can be obtained 
from the non-equilibrium force between two spheres [TH] , 
by appropriate summation in the optically dilute limit. 

We start with F 2 , where we get from Eqs. Q and 
^ in the limit of Ej — > 1, (note that this means that Sj 
is much larger than the radius of the objects) 
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On the other hand, Eq. (6) from Ref. [TH] for the in- 
teraction force between two spheres can be expanded in 
the dilute limit to yield the force between two volume 
elements V\ and U 2 , 2 



The force between two cylinders follows then as an inte- 
gral over the volume of the two cylinders, see Fig. [3j 
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Using this and Eq. (18), the terms in Eqs. (16) and (17) 
are recovered. Since the consistency of the two results is 
physically expected, it provides a useful check for both 
calculations. We note, however, that we could not recover 
the term proportional to d~ 2 in F^ 1 ' (corresponding to 
d~ 3 in Eq. (18)) due to mathematical difficulties in the 



integral over special cylindrical functions. 

We note also the corresponding equivalence for the self 
force, though somewhat more involved mathematically: 
In Eq. (7) from Ref. [15], the self force decays as dr 2 
with an oscillating prefactor exp (2iu)d/c). Pairwise sum- 
mation of this result yields exactly a d -3 / 2 separation 
dependence modulated by exp (2iujd/c), as discussed in 
relation to Eq. (15 1. Moreover, if we represent a plate of 



finite thickness by an infinite set of cylinders and place 
a small sphere in front of this plate, the self- force act- 
ing on the sphere scales as d^ 1 with the prefactor of 
exp (2iojd/c). This matches perfectly with the result for 
the sphere/plate geometry, in Eq. (19) from Ref. [TS], in 
the limit of large separations. Thus, in the dilute limit 
one regains the expected connections between different 
geometries. 



III. NUMERICAL EXAMPLES 

In this section we numerically compute the total 
Casimir force between two cylinders, in one case com- 
posed of dielectrics (SiC) and of conductors (tungsten) in 
the other. The equilibrium contributions to the force (see 
Eq. ([T])) were computed using the methods of Rec. [39] 
as implemented in the SCUFF-EM code suite [43] , 



Silicon carbide 



Figure 1(a) depicts the forces on SiC cylinders of radii 
R l= R 2 = 0.1/xm in both cold (0 K) and warm (300K) 
environments. The optical properties of SiC are modeled 
by the dielectric function [44], 
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2 We have an additional minus sign compared to Eq. (6) from 
Ref. 1 181 . and we exchanged indices 1 «-» 2. 



J LO 
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where = 6.7, ujlo = 0.12eV, ujto = 0.098eV, 
7 = 5.88 x KT 4 cV. We evaluated Eqs. ^ and (pj) 
numerically, restricting n and m to orders {—1,0,1}, 
and omitting quadratic terms in T, which is justified for 
{Rj} <C {At. }j {Sj}, where the forces per unit length will 
be proportional to R\R\. We note that Rj = 0.1/im is 
roughly an upper bound for the validity of this asymp- 
totic behavior, as we have checked by including higher 
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order terms. Due to the small thickness of the cylin- 
ders, the resulting phenomenology is in close analogy to 
the case of two spheres [T5] (compare also Sec. II D I: 
In the case when T env =0K, the force starts to devi- 
ate strongly from its equilibrium value around d ~ At/2, 
where At = ~ 7.6/xm. Cylinder 1 is repelled at large 
d if T 2 = 300K due to radiation pressure with a force 
that decays with distance as d^ 1 . On the other hand, if 
additionally T\ = 300K, the oscillating force is ap- 
preciable and in fact dominates the total force for large 
d if T\ — OK; the net force now has many zero crossings, 
where every second one is a stable point of zero force. 
The wavelength of the oscillations is roughly 6 (im due 
to the optical resonance of SiC at wavelength Ao ~ 12/im 
(the length of the optical path from cylinder 1 to cylin- 
der 2 and back is 2d, and for a sharp resonance of the 
dielectric function of cylinder 1 at Ao, we have construc- 
tive interference at 2d = Ao,2Ao, . . . , for an oscillation 
wavelength of Ao/2). Note that this figure also provides 
complete information about the force on cylinder 2: e.g., 
in case {Ti = 0,22 = 300K}, the blue curve shows the 
force acting on cylinder 1, while the red one represents 
the force on cylinder 2. At the crossing of the solid red 
and dashed blue curves the two cylinders feel equal forces 
in the same direction, an effect which might have less 
practical importance for tubes compared to spheres |18j . 



B. Tungsten 

Let us turn to conductors at high temperatures. Fig- 
ure [2] illustrates the total Casimir force for two very 
thin tungsten cylinders (nanowires) of equal radii Rj = 
0.02/im, and large temperature differences of 2400K. Wc 
have chosen this value for Rj as it corresponds to the 
maximum in emissivity of an isolated tungsten cylin- 
der [37] (which is connected to the skin depth being of the 
same order). As the non-equilibrium force strongly de- 
pends on the heat radiation of the objects, we expect it to 
be also comparatively large for this value. Also, tungsten 
has a relatively high melting temperature, and the value 
of 2400K promises large effects. Note that Eqs. and 
^ are not valid here, and the force is not proportional 
to R\R\, despite the small thickness of the wires. 

The following dielectric function for tungsten was 
used H5], 



2 



2ttcco ^ X rq 

q=l H 



(20) 



dielectric function of tungsten at T = 2400K H5] 3 . 

The relative deviation of the total Casimir force from 
its equilibrium form is not as pronounced as in Figjl] 
since equilibrium forces between conductors decay much 
slower with distance compared to dielectrics. In partic- 
ular, a Drude dielectric function gives rise to an equilib- 
rium force scaling as d~ A / log (d/R) [6], in good agree- 
ment with our numerical results. (For dielectrics, our 
numerical result suggests a d~ 7 law of the equilibrium 
Casimir force in the small separation retarded regime). 
Nevertheless, as in the dielectric case, the force in Figj2] 
starts to deviate from the equilibrium curve at approx- 
imately At/2 with a corresponding At ~ 0.95^m. In a 
cold environment, T env = OK, the force is attractive at 
small separations, whereas at large distances it becomes 
repulsive if Ti = 2400K. When we increase the environ- 
ment temperature to T env = 2400K, the force shows an- 
other feature for T\ = T 2 — OK, where it is attractive 
for both small and large distances, but for an interme- 
diate region it becomes repulsive yielding a stable point 
of zero force in the vicinity of d ~ 4/xm. This intermedi- 
ate repulsion is due to the contribution from evanescent 
waves (which is repulsive because i* 1 - enters Eq. ([!]) with 
a minus sign), and then for higher d propagating waves 
make the force attractive again. In contrast to the di- 
electric displayed in Fig. 1(b) [ no oscillations are visible 
in Fig. [2] We attribute this to the absence of sharp peaks 
in the dielectric function of tungsten. 

The non-equilibrium force in Fig.[2]is much larger com- 
pared to Fig.[TJ again, as a consequence of both replacing 
the dielectric materials with conductors, as well as choos- 
ing the higher temperature T — 2400K instead of T — 
300K. At the point where the total Casimir force deviates 
from the equilibrium force the magnitude of the force is 
roughly 10 6 times larger for tungsten compared to SiC- 
approximately 5 x 10 6 larger per unit mass. For tung- 
sten nanowires at T = 300K (we do not provide the plot 
for this case) the non-equilibrium effects are still three 
orders of magnitude larger when compared to SiC. Thus, 
the strong enhancement is a result of both material prop- 
erties and high temperatures. Furthermore, considering 
wires of length 1/mi at high temperature T = 2400K, the 
force at the point when non-equilibrium forces start to 
dominate is approximately lOfN (see Fig. [2), significantly 
larger than the weight F G ps 0.24fN (also depicted in the 
figure) of the nanowires. 

There is another force we can use for comparison: If 
an electric current is used to heat the wire [14] [46] , one 
requires a current of 17/xA to maintain the temperature 
T = 2400K for tungsten nanowires (obtained by equating 



where A is the wavelength in vacuum, c is the speed 
of light and eo is the permittivity of vacuum in SI 
units. The remaining parameters are : o\ — 1.19 x 10 6 
ohm m j er 2 — 0-25 x 10 6 ohm m , A rl = 3.66/im, 
A r 2 = 0.36/xm. These values were obtained by fitting the 



3 For the accurate evaluation of the total Casimir force with 
Eq. JH , one must use the cylinder dielectric function at its cor- 
responding temperature. Thus, Eq. (|20| is strictly only appro- 
priate for Ti = T2 = 2400K in Fig. |2] However, the dielectric 
function does not depend significantly on T and the curves for 
Tj = OK in Fig. [2] are good estimates. 



heat losses described by Joule's law to the heat radiation 
for a signle cylinder predicted in Ref. [37]). Such cur- 
rents correspond to an Ampere's force of about 0.15fN 
for d = 0.4/mi, which is two orders of magnitude smaller 
than the nonequilibrium correction to the Casimir force 4 . 
Thus, when switching on an AC or DC current in the 
wires, the change in force between them due to heating- 
induced non-equilibrium Casimir force is much larger 
than that due to Ampere's force. We thus conclude that 
hot conducting nanowires are promising candidates for 
measuring or using non-equilibrium Casimir effects. 



Appendix B: Small R expansion of the T operator of 
the cylinder 



In order to derive Eqs. 0, (Til and (III) we need 
the expansion of the T operator in terms of ujR/c. For 
a cylinder made of isotropic material with magnetic per- 
meability and dielectric permittivity e{uj), we find 
for the limit R < {5, c/uj} [37], 
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Appendix A: Cylindrical harmonics 

Following Ref. [47], the EM cylindrical harmonics can 
be written as, 



rpMM _ rpM 

1 l,k z " -1 



mm _ «r kl( P - l)(e + l) + (fi+ 1)(b - 1) 
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x (ujR/c) 2 , 
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where H n is the Hankcl function of the first kind of order 
n. M B) t z and ~N n> k z correspond to outgoing magnetic 
multipole (TE) and electric multipole (TM) waves re- 
spectively. Also, k z and q are the wavevectors parallel 
and perpendicular to the cylindrical z-axis respectively, 
satisfying the relation q = y?k 2 — k z , k — u>/c. Hn 1 ^ cor- 
responds to the first derivative with respect to the argu- 
ment. Furthermore, we denote the corresponding regular 
waves by RM n ^, and RN^n,, which differ from regular 

ones by replacing Hn with the Bessel function J n . 

The above solutions correspond to transverse waves, 
i.e. V-M„ ifcz = V-N nifez =0. Moreover, they obey useful 
relations M n ,^ = fjV x N„, ls) N„, fc> = x M n ,k z - 
These relations are also valid for outgoing waves. 



4 This estimate neglects possible couplings between AC or DC cur- 
rents and thermal current fluctuations |48l 



rpMN _ rpNM _ rpMN 

1 l.k z ~ 1 l,k z — 1 -l,k, 



where k z — k z /k. 
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2(e+l)( M +l) 



x (ljR/c) 2 , 



(B.5) 



Appendix C: Auxiliary functions 

In defining functions below, the superscript emphasizes 
that we deal with the interaction force, whereas the nu- 
merical subscript indicates the force's power of decay in 
the axis-to-axis separation between cylinders. 



l (£iM,£ 2 H) 



45 
2048 



Im 



1 



e 2 + 1 



1 



x [(| £l | 2 - 1)(4(33 + 5 Re[ei]) + (7 + 3Rc[e 1 ])|e 2 + 1| 2 ) 
+(Re[ £l ] 2 -l)(40 + 6|e 2 + l| 2 )] ; 

(C.l) 



5r(£l ' £2)= 256 Im 



1 

e 2 + 1 



1 

ki + ii 



x [Im[e 1 ] 2 (|e 2 + 1| 2 (7 - Re[ £l ]) + 12Rc[ £l ] + 76) 

+(Re[ei] 2 - l)(|e 2 + 1| 2 (5 - Re[ei]) + 12Re[ei] + 100)] ; 

(C.2) 
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g i 1 n (e 1 (u),e 2 (u)) = — Im 
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Im 
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(|ei + 1| 2 N + 1| 2 + |ei + 1| 2 + |e 2 + 1| 2 + 36) ; 
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